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Determination of airway resistance and lung compliance 



(57) A method and apparatus is disclosed for deter- 
mining airway resistance and lung compliance using an 
electrical circuit model (36) wherein at least one com- 
ponent parameter is non-linear. The system (10) non- 
intrusively obtains pressure and flow data signals 
(20,24) from a pressure transducer (22) and a laminar 
flow element (18) without interrupting or interfering with 
normal breathing and gas supply to a patient (12). An 
invariant exponential is determined empirically based 



on physical characteristics (14) of the airway The non- 
linear airway resistance and lung compliance can then 
be calculated based on the sensed flow rate, gas pres- 
sure, a calculated gas volume, and the invariant expo- 
nential using linear techniques. The resulting airway re- 
sistance can be normalized to a standard reporting flow 
rate value. The system is particularly useful in anesthe- 
sia applications, but is also useful in any breathing sys- 
tem where fresh gas is supplied constantly from a gas 
source other than a ventilator. 
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Description 

BACKGROUND OF THE INVENTION 

[0001] The present invention relates generally to determining airway resistance and lung compliance, and more 
particularly to using a circuit model approach to find the ainway resistance and lung compliance where at least one 

component is non-linear, . 

r0002] In an anesthesia procedure, it is advantageous to know with some certainty the airway resistance and lung 
compliance in order to ascertain the suitability of the ventilatory management system. Although many attempts have 
been made to determine airway resistance, most assume the resistance to be linear, when in fact, it is not. Others do 
not consider the effects of lung compliance. Therefore, both types of systems fall far short of determining either pa- 
rameter with any certainty and in some cases, results in gross inaccuracies. 

r0003] It will be shown herein that those systems that assume a linear resistance relationship between pressure and 
flow in the patient ainway vwll not function properly on intubulated patients. It is also believed that since an endotrachea^ 
tube follows the natural path and shape of a patient's ainway. such linear systems will not function properly if applied 
directly to the patient's airway passage. The error found in the results from such linear techniques for determining 
resistance and compliance will increase dramatically virith varying ventilation flow rates. Changing flow rates is common 
in anesthesia procedures and can be caused by changes to ventilatory settings of tidal volumes, inspiratory flows, or 
fresh gas delivery to the breathing circuit. In order to assume such a linear relationship then, one must maintain a 
constant flow rate. However, such an undesirable dependence on requiring a constant flow rate, can also impose errors 
to the very parameters being estimated because flow rate changes during each breath cycle as well. 
[0004] Some prior art systems require the injection of an excitation flow into the breathing circuit or an inspiratory 
pause in order to calculate the ainway resistance and lung compliance. However, such techniques are not practical 
during an anesthesia procedure. In fact, in any breathing system where fresh gas is supplied constantly from a gas 
2S source other than the ventilator, an inspiratory pause cannot be imposed. 

[OOOS] Other known systems use a forced high frequency oscillation to determine airway resistance and lung com- 
pliance The problem with this system is that no one good resonant frequency can be determined for all patients. 
Attempting to find the correct frequency for each patient would be time consuming and not practical, 
[00061 One early attempt at determining lung airway resistance non-linearly is disclosed in U.S, Patent 3,036.569. 
However merely finding a resistance at one flow rate does not provide sufficient data to be useful in anesthesia pro- 
cedures 'it has also been found that pressure is not a function of resistance only, but also of compliance. Further, this 
reference requires an apparatus that forces air into the lung in order to perform the calculation, which would interfere 
with normal breathing and with anesthesia flow. . . ^ , 

[00071 It would therefore be desirable to have a system, including a method and apparatus, that does not interfere 
with normal breathing, is non-intrusive with the normal flow and pressure during an anesthesia procedure, can measure 
pressure and flow on expiration, does not interrupt or interfere in any way with the respiratory pattern, and can find 
both ainway resistance and lung compliance, while still reporting the non-linear air resistance at a standardized flow rate. 
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SUMMARY OF THE INVENTION 



[0008] The present invention provides a system for determining the non-linear ainway resistance and lung compliance 
using a circuit model approach that overcomes the aforementioned problems. 
[00091 In accordance with one aspect of the invention, a non-linear method of establishing airway resistance and 
lung compliance is disclosed using a circuit model. The method includes the steps o1 sensing a gas flow rate through 
45 an ainway and sensing a gas pressure in the airway, then calculating a gas volume from the gas flow rate, and deter- 
mining an invariant exponential based on physical characteristics of the ainway Airway resistance and lung compliance 
can then be accurately calculated based on the gas flow rate, the gas pressure, the gas volume, and the invariant 
exponential at any flow rate. . 
[00101 in accordance with another aspect of the invention, an apparatus is disclosed to determine airway resistance 
so and lung compliance. The apparatus includes an airway capable of communicating external gas with a patient's lungs, 
a gas flow rate sensor attached to the ainway to sense a gas flow through the airway and produce a flow signal in 
response and a gas pressure sensor located in the airway to sense a gas pressure across the airway and produce a 
pressure signal in response. A processor, such as a computer, a central processing unit, a microcontroller, or any other 
type of processing unit, is connected to the gas flow and pressure sensors to receive the flow and pressure signals. 
55 The processor is programmed to calculate airway resistance and lung compliance using a non-linear model having at 
least one non-linear component. 

[0011] The current system does not require injection of excitation flows into the breathing circuit or an inspiratory 
pause to calculate the airway resistance and lung compliance, as required in the prior art. The present invention ac- 
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quires pressure and gas flow rate signals from pressure transducers which measure the relative airway pressure and 
the pressure across a laminar flow element without interruption or interference with normal gas flow 
[0012] The analog signals acquired are digitized and supplied to a processor which is programmed to calculate the 
non-linear airway resistance and lung compliance according to an electrical circuit model. The total airway pressure 
is equal to the sum of the pressure due to the flow rate, the pressure due to the volume, a pressure constant, and the 
pressure due to flow change, which in most cases is minimal and can be ignored. The pressure due to flow rate, or 
resistance, Is empirically determined to be non-tinear. Various airways such as endotracheal tubes are analyzed and 
are fit to a common equation in order to determine an invariant exponential. Data is acquired at three convenient points 
along the expiration cycle in order to solve the circuit model equation by common linear algebra techniques. The 
resistance can then be converted to a standardized reporting flow rate. Such conversion however is optional since the 
resistance calculated is accurate at any flow rate, unlike linear resistance models. 
[0013] Various other features, objects, and advantages of the present invention will be made apparent from the 
following detailed description and the drawings. 

15 BRIEF DESCRIPTION OF THE DRAWING 

[0014] The drawings illustrate the best mode presently contemplated tor carrying out the invention. 

[0015] In the drawings: 

[0016] Fig. 1 shows a block diagram of a system in accordance with the present invention as applied to a human 
20 subject. 

[0017] Fig. 2 is a circuit schematic of the circuit modeling approach encompassed in the system of Fig. 1 . 

[0018] Fig. 3 is a graph showing flow versus pressure across the resistor for various airway tubes. 

[0019] Fig. 4 is a flow chart used in implementing the system of Fig. 1 . 

25 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

[0020] Fig. 1 shows a system 10 which includes an apparatus to determine ainway resistance and lung compliance 
in a patient 12. In one application of the present invention, an endotracheal tube 14 acts as an airway between an 
external oxygen source (not shown), which can include an anesthesia component, and the lungs 16 of patient 12. The 

30 system 10 includes a gas flow rate sensor 18 attached to the airway 14 to sense a gas flow therethrough and produce 
a flow signal 20. A gas pressure sensor 22 is also located in airway 14 to sense a gas pressure therein and produce 
a pressure signal 24 from the sensed gas pressure in airway 14. The flow and pressure signals 20, 24 may be provided 
by pressure transducers in the flow sensor 18 and the pressure sensor 22 which measure the relative pressure across 
a laminar flow element and the relative airway pressure, respectfully. The signals 20, 24 are proportional to the pressure 

35 and flow and are filtered to remove noise and errant signals by the A/D and signal conditioner 26, which also converts 
the analog signals to digital form for processing by a CPU 28. The data acquisition occurs on a discrete time basis; 
that is, the A/D converter 26 establishes a data value for the respective signal over its sampling interval, later referred 
to as the sampling time. 

[0021] The CPU 28 is connected to the gas flow sensor 18 and the pressure sensor 22 via the A/D converter 26 to 
40 receive digitized flow and pressure signals 30. The CPU 28 is programmed to calculate the ainway resistance Rp and 
the lung compliance Cl using a non-linear circuit model having at least one non-linear component, as will be further 
described with reference to Figs. 2-4. Once the airway resistance Rp and the lung compliance Cl are known, repre- 
sentative signals 32 can be transmitted to an external monitoring apparatus 34 to monitor the ventilatory management 
system. Although Fig. 1 shows the system 1 0 ot the present invention applied to an airway in the form of an endotracheal 
45 tube 14, the present invention is not so limited. For example, the patient may be fitted with a mask having an external 
airway 14 attached thereto. 

[0022] Referring to Fig. 2. a circuit schematic 36 is used to illustrate the circuit model approach used to determine 
the airway resistance Rp and the lung compliance Ci_ according to the present invention. To more accurately find the 
airway resistance Rp and the lung compliance Cl, the mechanics of the respiratory system are approximated by a 
50 mathematical model that relates airway pressure P^w at the vicinity of the patient's mouth to the bidirectional gas flow 
in and out of the lungs. This mathematical model is illustrated as the circuit model of Fig. 2. The functional relationship 
is described as: 
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^aw + + P"" (volume) + P2{flow) + P3(flow changes) (1) 
[0023] Each of the Px terms on the right side of Equation (1 ) are various pressure contributions from the pulmonary 
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T^Z^ In T^^^^ the functional relationship of the restriction term P2, indicated by resistor 44^ .s 

5 Omr^ ?ube cTve 54 shows that relationship for a 6.0mm. tube, and curve 56 shows the non-linear properties for a 
S an exponential function, as shown in Equation 2, or a polynomial senes. as shown in Equation 3. 



P2(F) = K1 M (t)" 



P2 (F) = AO = A1*f(t) + A2*f (t) 



2 + ... 



30 

r002G] in a typical ventilatory range during anesthesia, the parameters Kt and '^^.^'^''T^'^'^;^^^^^^ 

aTv sinqle breath The f (t) term is the instantaneous bidirectional flow. From the graph .n Fig. 3 the data for t^^^ pressure 

Td Z can then be used to find the unknown exponential n. For the various endotracheal tubes shown in F.g. 3. 

the following equation is used to model and calculate patient pulmonary mechanics: 

P3, = L.1/C,*v(t).KpM(tr W 

[0028] Equation 4 is the particular equation that relates to the general Equation 1 with '^fJl'^'l'i^Zo^^^^^ 
P3 assumed to be zero The L term is the PO constant term. C, is modeled as the compliance o he lung and v t) 'S 
- thetnlaranlous volume in the lung. The lung volume v(t) is found by integrating the bi^^^^^^^^^^^^ b 
described in more detail with reference to Fig. 4. The product of the inverse lung ^^^^^""^f/"'!'^^^;^^^^^^^^ 

?M29r T^e curves ot F°g 3 were plotted by placing a constant gas flow through each tube and measuring the flow 
riput as w^as he change in pressure across the tube. The invariant exponential is found by fitting each curve 
to a crmo^f ^cti^n and although t'he K, term changes ,or each tube, or resistor, the 2;P0--^ J-^^Vc! ^.^^ 
in this case and it is presumed tor all tubular airways, the invariant exponential is 1.7 which ^^P'®^^"ts me j^^^^^^ 
n mL endoUacheal tube. It is understood that different geometries of airway configurations may '^^^"f '"^^^^^ 
exponenSal However, during anesthesia, the endotracheal tube dominates the ain.ay resistance. These tubes are 
Snar to The those used empirically to derive the exponent. In masked cases where patients are not intubated, the 
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endotracheal remains to be the dominant airway resistance. Consequently, the invariant exponent value of 1 .7 applies 

l^To? trnSw be'evldent. having a value for the exponent, and measurements tor the bidirectional flow rate^ f 
m Z airway pressure, P^^. and a calculated volume. v(t), the calculation of the airway res.stance, Rp, and the lung 
lompnaTce C, is reducedTo a problem of value identification for C, and K^, and ultimately the ^^^^^[^^^^"^^ ^^^^^^ 
R in the p eferred embodiment, the solution presented uses simultaneous equat.ons of three sets of data po,n s to 
Slve for th'e unknown C, and K, Specifically, three convenient points are chosen obt^'" ^^^^ ^he f,rs^^^^ ^^'^ 
T when the flow rate is equal to zero at the beginning of an expiration. The second is at a time when the flow rate 
•s^aTa minimum afir time\, , and the third is at a time T3 after T, when the flow rate °' ^^^"^.^^^^^^ 

words, the three sets of data points are taken at the end of inspiration, at minimum expiratory flow and at 50./o expiratory 

l^'l] in this case, expiration, or flow out of the patient, is chosen negative. The equations ^"^^^^f f g"^^^^^^^ 
hen be represented in matrix notatbn and may be solved by various known techniques. For example a basic matro^ 
auqmentXn and row reduction approach can be used for simplicity. However, one skilled ,n the art will recognize that 
: 'ous Ssr t^hniques can be ir^p.emented to solve for the unknown lung compliance ^^^^^^^^^^^^^^^ 
resistance Kp. such as regression or digital filtering. Such methods are less sensitive to measurement noises but are 

?3TpS"cruse?s are familiar and comfortable with a resistance representation, that Nnearly relates airway 
p^ssure and flow rate. To meet this expectation, all the non-linear ainA,ay resistances are mapped to hnear resistances 
rlferenced a, a standardized gas flow rate before it is reported. This linear airway resistance vanes with flow rate and 
Lhol on?y be compared at the referenced flow rate. The following relationship is used to report the airway resistance 

at a relerenced flow rate: 
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n-1 (5) 



where Rp is the linear airway resistance at a referenced flow rate f,,. In practice, it is convenient to report the airway 

resistance at a standardized 30 liters per minute flow rate. 

S Although the preferred embodiment describes the aforementioned relationships for an anesthesia app heat on. 
irprLem inveSion is' readily applicable to other ventilatory conditions or environments wherein the terms contributing 
to ainway pressure can be described by different relationships or constants. includes 
r0O341 Referring to Fig. 4, the software algorithm is described in flow chart fomn. The flow chart of Fig. 4 includes 
Sa acquisL at three'points during expiration, volume determination, calculation of the unknowns, -d conversion 
of the ai^rway resistance to a standardized flow rate. Upon start up 58 all values are initialized to 1 60 and 62. The 
anlloq values for the pressure and flow are read 64 from the pressure and flow sensors and the ana og signals are 
?he?dlitted 66 The Sow and pressure values at the minimum flow (F_ZERO) are detem^ined by cominuously mon- 
tor^g tt present and previous'flow rates to differentiate between inspiration and expiration 68. 70. When the presen 
ow (F NEW) is zero or less than zero, and the previous flow (F.LAST) is above zero, then the minimum flow (F MINI) 
S ie-en found 72. indicating the beginning of an expiration cycle in which the values for flow and Pressure can ^ 
determined and saved as the minimum flow values (V_ZERO. P_ZERO) 74. If the present flow (F NEW) 's at zero 
fhen the £w volume, and pressure values are simply saved. However, if the Present flow is less t^-zera then the 
values are interpolated for zero flow and the interpolated values are saved for V_ZERO and P_ZERO. Once in an 
expiSircyde and zero flow has not yet been reached, the volume value is updated 76 by adding the previous value 
for the volume to the product of the latest flow value and its respective sampling time. .„n,in,. 
S)351 Next a minimum flow determination is made. After the zero flow value has been found, the system cont.nu- 
oT^v monies fhT flow signal to determine when it has reached a minimum value 78. This is accomphshed by cont.n- 
uSy corJparing t^^ present flow value (F_NEW) with the previously saved minimum (F_MIN). When the present 
^alue is less than the previous minimum, then the minimum is set to this present value at 80 and the volume and 
p^ssureTor tSlow ra^ value are saved as V_MI N and P.MIN. Again, in determining these values il has earlier been 
assumed that the flow rate out of the patient is negative. The flow rate could be assumed positive, with corresponding 

r ird^^r rd^e^termined at a SO. flow rate. To find F^SO. the system ^^^^^^^^ ^ 
low signal to determine when it reaches 50% of the previously found minimum .^^-/^^^^j^^^ 
compllshed by comparing the present flow value (F_NEW) wrth 50% of the minimum « f -^^^^ f^^^^hen 

55 the present value is less than half the minimum value 84. the volume and pressure related to this 50 /» flow rate value 

^SaT'l^nrend ora^expSion cycle 88. 90, the unknowns Kp, C, and L can be found at 92, as previously set 
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forth. The resistance is then standandized 94 and can be reported to an external monitoring apparatus 96 and the 

rST TXlT^se!L, be iost fror. the lung thereby making the lung volume actually smaller than the ^e- 
qTd bidreSal low. The to'al volume loss within a breath can be determined by the difference of the 'nsp-red t.da 
vSime 10 he expired tidal volume. The Instantaneous volume losses may be estimated by apportioning the ratio of 
thttma I vcTume loss In that breath to the instance of volume measurement. The ratio would "^^^^f ^"^'"^^^^'3"^;'^ 
This would then minimize the affect of volume loss in the resistance compliance calculation o1 the present invent.on^ 
■2^397 Accordingly, the present invention also Includes a non-linear method of establishing ^jrway resistance and 
u^q comphan e using a circuit model. The method includes the steps of sensing gas flow rate through an airway and 
sensing gTpressure in the airway. The method also Includes calculating a gas volume ^om the gas f'ow rate a d 
dete^n?ig an Invariant exponential based on the physical characteristics of the ainA,ay. A-ay res.tance and K^ng 
compliance can then be calculated based on the gas flow rate, the gas pressure, the gas volume, and the invariant 

AS ^Z^Z:%;:Ze to Fig. 4, the step of calculating gas volume Includes —^^^^ 
Lnd inspiration flow rates and multiplying each sensed expiration gas flow rate by a corresponding san^pling time for 
Tcu S qas voLe sample. The results are then integrated as a series of current gas volume samp es during the 
expirron'cyc'e After at least three sets of data are acquired, the airway resistance and lung corripUartce can be 
IZutZ by either forming a matrix of the acquired data and soK/ing the matrix, or with the use of regressive techniques 

TeTesem InTention has been described In terms of the preferred embodiment, and it ■s/ecognized that 
L^JlJLnts. attematives, and modifications, aside from those expressly stated, are possible and within the scope of 
the appending claims. 



25 Claims 

1. A nonlinear method o1 establishing airway resistance and lung compliance using a circuit model comprising the 

Steps ot: 

30 sensing gas flow rate through an airway; 

sensing gas pressure in the airway; 
calculating a gas volume from the gas flow rate; 

determining an invariant exponential based on physical characteristics of the ainway ; and 
calculating airway resistance and lung compliance base on the gas flow rate, the gas pressure, the gas volume 
35 and the invariant exponential. 

2. The method of claim 1 wherein the step of sensing gas flow rate and pressure are further defined to occur in an 

endotracheal tube. 



40 



3. The method of claim 1 wherein the step of calculating gas volume further compnses the steps ot: 

differentiating expiration and inspiration flow rates; uoinmP 
multiplying each sensed expiration gas flow rate by a corresponding sampling time for a current gas volume 

sample; and 

45 integrating a series of current gas volume samples during an expiration cycle. 

4. The method of claim 1 further comprising the steps of: 



so 



55 



acquiring at least three sets o1 data trom the sensing steps; and 

forming a matrix of the acquired data for calculating the airway resistance and lung compliance. 

5. The method of claim 1 further comprising the steps of: 

acquiring at least three sets o1 data from the sensing steps; and 

regressively calculating the airway resistance and lung compliance using the acquired data. 

6. The method of claim 1 wherein the airway is tubular and the invariant exponential is approximately 1 .7. 
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7. The method of claim 1 wherein the step ot calculating airway resistance and lung compliance is lurther defined as 

solving the equation: 

P^^ = L-H/CL*v(t) + KpM(t)" 

where P is the sensed airway pressure, v(t) is the calculated airway volume, f (t) is the flow rate, n is the 
invariant exponential, L is a constant term, C, is the lung compliance, and K, >s the airway resistance. 

8. The method of claim 7 wherein the ain«,ay resistance is normalized to a standard flow rate given by; 



p p 



where F.e,"- is a referenced flow rate and Rp is the normalized airway resistance. 

9. The method of claim 1 wherein the step ot determining an invariant exponential is dependant on a shape and size 
o1 the airway and is determined empirically. 

10. An apparatus to determine airway resistance and lung compliance comprising: 

r^rwTeri^^^^^^^ riTrtSrethrough and produce a flow signal 

TgaJ pre;sure sensor located in the airway to sense a gas pressure therein and produce a pressure signal 

rrocrsor connected to the gas flow and pressure sensors to receive the flow and pressure signals, the 
pr^essor "rogr^^^ to calculate airway resistance and lung compliance using a non-hnear model having 
at least one non-linear component. 

11 The apparatus of cteim 1 0 wherein the non-linear model is based on a circuit model approach and a characteristic 
' of resistance includes an invariant exponential in the non-linear component, 

1 2 The apparatus of claim 1 0 wherein the processor is programmed to accept an external invariant exponential input, 
the invariant exponential based on a particular ainway configuration. 

1 3. The apparatus of claim 1 1 wherein the airway is tubular shaped and the invariant exponential is approximately 1 .7. 

14. The apparatus of claim 10 wherein the processor is further programmed to solve the equation: 

P^=L + 1/CL*v{t) + Kp*f(t)" 

Where P is the sensed airway pressure, v(t) is the calculated airway volume, f(t) is the flow rate, n is the 
invariant exponential, L is a constant term. C, is the lung compliance, and is the a,n.ay resistance. 

15. The apparatus of claim 14 wherein the airway resistance is normalized to a standard flow rale given by: 



ref 



Where f,^,^-' is a referenced flow rate and Rp is the normalized airway resistance. 
16. The apparatus ot claim 10 further comprising: 

a signal conditioner to filter errant and noise signals received from the sensors; and 

an A/D converter connected to receive analog signals from the gas flow rate sensor and the gas pressure 
sensor and produce digital signals in response. 
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'l7. The apparatus of claim 10 wherein the processor initially calculates a gas volunae from the received signals. 

18. The apparatus of claim 17 wherein the processor calculates the gas volume by. 

=i differentiatina expiration and inspiration flow rates; 

muSng each sensed expiration gas flow rate by a corresponding sampling time for a current gas volume 

sample; and . 
integrating a series of current gas volume samples during an expiration cycle. 

10 19. The apparatus ot claim 10 wherein the airway is an endotracheal tube. 
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FIG. 4 
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